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ABSTRACT 



Context. The plateau in the abundance of 7 Li in metal-poor stars was initially interpreted as an observational indicator of the pri- 
mordial lithium abundance. However, this observational value is in disagreement with that deduced from calculations of Big Bang 
nucleosynthesis (BBN), when using the Wilkinson microwave anisotropy probe (WMAP) baryon density measurements. One of the 
most important factors in determining the stellar lithium abundance is the effective temperature. In a previous study by the authors, 
new effective temperatures (T^) for sixteen metal-poor halo dwarfs were derived using a local thermodynamic equilibrium (LTE) 
description of the formation of Fe lines. This new T B g scale reinforced the discrepancy. 

Aims. For six of the stars from our previous study we calculate revised temperatures using a non-local thermodynamic equilibrium 
(NLTE) approach. These are then used to derive a new mean primordial lithium abundance in an attempt to solve the lithium discrep- 
ancy. 

Methods. Using the code MULTI we calculate NLTE corrections to the LTE abundances for the Fei lines measured in the six stars, 
and determine new 7Ws. We keep other physical parameters, i.e. log g, [Fe/H] and constant at the values calculated in Paper I. 
With the revised r c(f scale we derive new Li abundances. We compare the NLTE values of r c(f with the photometric temperatures of 
Ryan et al. (1999, ApJ, 523, 654), the infrared flux method (IRFM) temperatures of Melendez & Ramirez (2004, ApJ, 615, 33), and 
the Balmer line wing temperatures of Asplund et al. (2006, ApJ, 644, 229). 

Results. We find that our temperatures are hotter than both the Ryan et al. and Asplund et al. temperatures by typically ~ 110 K - 
160 K, but are still cooler than the temperatures of Melendez & Ramirez by typically ~ 190 K. The temperatures imply a primordial 
Li abundance of 2.19 dex or 2.21 dex, depending on the magnitude of collisions with hydrogen in the calculations, still well below 
the value of 2.72 dex inferred from WMAP + BBN. We discuss the effects of collisions on trends of 7 Li abundances with [Fe/H] and 
T c ff , as well as the NLTE effects on the determination of log g through ionization equilibrium, which imply a collisional scaling factor 
Sh > 1 for collisions between Fe and H atoms. 
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1. Introduction 



Since its discovery bv lSpite & Spitd ([T982), many studies of the 
plate au in lithium in metal-poor dwarfs have been undertaken, 
e.g.lSpite etal](ll996l).lRvan et al.ld2000h,lMelendez & Ramirez! 
(2004), Bonifa cio et all d2007l) and lAoki et alj d2009l) . confirm- 
ing its existence. Most studies find a comparable Li abundance 
(A(Li£] w 2.0 - 2.1 dex) yet discrepancies still exist, in particular 
the high value found by iMelendez & Rami rez (2004) (A(Li) = 
2.37 dex). However, the biggest discrepancy comes from a com- 
parison of the primordial abundances inferred from observations 
and that derived from Big Bang Nucleosynthesis (BBN) with the 
WMAP constraint on the baryon density fraction, Q.B,h 2 , which 
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leads to A(Li) = 2.72 dex dCvburt et alj|2008l) . This is what has 
become known as the "lithium problem". 

Several possibilities have been proposed to explain this dis- 
crepancy. Broadly these are: systematic errors in the derived stel- 
lar Li abundances; errors in the BBN calculations due to un- 
certainties in some of the relevant nuclear reaction rates; the 
destruction of some of the BBN-produced Li prior to the for- 
mation of the stars we have o bserved; the introduction of new 
physics that may affect BBN dJedamzik & Pospelovll2009l) ; or 
the removal of Li from the phot ospheres of the star s through 
their lifetimes (see introduction to lHosford et al.ll2009L Paper I, 
for more details). The possible explanation under study in this 
work is that of systematic errors in the effective temperature 
(Jeff) scale for metal-poor stars. The effective temperature is the 
most important atmospheric parameter affecting the determina- 
tion of Li abundances. This is due to the high sensitivity of A(Li) 
to Jeff, with dA/dT e ff~ 0.065 dex per 100 K. One reason for the 
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spread in the observed A(Li) is the diffe rences in th e r e g s cales 
used by different au thors. For instance, ISpite et a"D (1 19961) and 
lAsplund et ail d2006l) derive a r eff of 5540 K and 5753 K for the 
star HD 140283, respectively. The scale of Mele ndez & Ra mirez 
(2004) is on average ~ 200 K hotter than other works. This goes 
some way to explaining their higher A(Li); other factors, such 
as the model atmospheres with convective overshooting used in 
their work, may also contribute to the discrepancy. It is impor- 
tant to confirm, or rule out, whether systematic errors in T e g are 
the cause of the Li problem, and in doing so address the need for 

other possible explan ations. 

In previous work dHosford et alj |2009. - Paper I), we utilised 
the exponential sensitivity in the Boltzmann distribution to x/T, 
where x 1S the excitation energy of the lower level of a transition. 
Using this, we determined r e ff's for eighteen metal-poor stars 
close to the main-sequence turnoff. This was done by nulling 
the dependence of A(Fe) on x f° r approx 80 - 150 Fei lines. 
Two T e ff scales were generated due to uncertainty in the evo- 
lutionary state of some of the stars under study. It was found 
that our temperatures were in good a greement with those de- 
rived by a Balmer line wing method bv lAsplund et al.l ( 2006) and 
those derived by photometric techniques bv lRvan et al.l (11999). 
However, our T e g scale was on average ~ 250 K cooler than 
temperatu res from the infrared flux me thod (IRFM) as imple- 
mented by iMelendez & Ramirez! (12004). This is not the case for 
all work done u s ing th e IRFM, the IRFM effective temperatures 
of lAlonso et al. I (fl996l) are similar to ours, for stars we have in 
common. 

The derived mean abundances in Paper I were A(Li) = 2.16 
dex assuming main-sequence (MS) membership and A(Li) = 
2.10 dex assuming sub-giant branch (SGB) membership. For the 
five stars that have a known evolutionary state, we calculated a 
mean A(Li) = 2.18 dex. It is clear that these values are not high 
enough to solve the lithium problem. However, the analysis of 
iHosford et alJ (2009) assumed that the spectrum was formed in 
local thermodynamic equilibrium (LTE). This is a standard way 
of calculating spectra, but oversimpli fies the radiat i ve tra nsfer 
problem, and it was acknowledged in Hosf ord et alJ (120091) that 
LTE simplification affect those results. Consequently, although it 
was shown that, within the LTE framework, systematic errors in 
the T e g scale are not the cause of the disparity between spectro- 
scopic and BBN+WMAP values for the primordial Li, we also 
need to assess the impact of non-local thermodynamic equilib- 
rium (NLTE) on the determination of stars effective tempera- 
tures. That is the aim of the current work. 

This work is not intended to be a full dissection of the meth- 
ods of NLTE, but rather an application of those more complex 
(and possibly more accurate) methods to derive a new T e s scale 
and to assess their impact on the lithium problem. However, to 
do this we need to delve, with some depth, into the processes 
of NLTE line formation, which we do in Sect. 2. This will give 
some understanding of the complexities and uncertainties that 
are involved and give the opportunity to make some generalisa- 
tions on the important aspects that need to be addressed. In Sect. 
3 - 5 we detail our calculations and results, and discuss these 
further in Sect. 6. 

2. NLTE Framework 

2.1. The necessity for NLTE 

With the availability of high quality spectra, the problem of cal- 
culating accurate chemical abundances often comes down to a 
better understanding of the line formation process. This is of 



particular importance to this work as the calculation of accu- 
rate level populations of the Fe i atom and source functions at 
the wavelengths of the Fe transitions is crucial to determin- 
ing r e (j from lines of different x- In LTE calculations, the level 
populations follow the Boltzmann and Saha distributions. These 
assume that the levels are populated, or depopulated, by colli- 
sional and/or radiative processes, that are characterised by the 
local kinetic temperature. In the deep layers of the atmosphere, 
at T5000 > 1, where T5000 is the optical depth at 5000A, LTE 
is a reasonable assumption. However, it tends to break down at 
optical depths T5000 < 1> i- e - through most of the line forming 
region of the photosphere. Therefore neglecting deviations of 
the level populations from LTE could lead to errors in the T^ff 
derived by excitation dependence. Furthermore, in NLTE calcu- 
lations, it is not only the level populations that differ from the 
LTE case. The radiative transitions of the atom must be explic- 
itly considered. The fact that the radiation field is no longer de- 
scribed by a Planck function, and certainly not a Planck function 
calculated for the local temperature, results in further changes of 
the spectrum relative to the LTE case. This last effect is very im- 
portant in metal-poor stars, where the reduced opacity/increased 
transparency of the atmosphere exposes shallow, cooler layers 
to the UV-rich spect rum coming from the deeper, hotter layers 
(lAsplund et al.lll999l) . 

For Fe in particular, different studies have come to differ- 
ent conclusions as to th e magnitude of the NLTE corrections. 
iThevenin & Idiartl d 19991) found that there can be corrections 
of up to 0.35 dex on Fei abundances for main-sequence stars 
at [Fe/H] as -3, and suggest that all work done on metal-poor 
stars s hould be carried out using NLTE methods. iGratton et al.l 
dT999l however, find negligible corrections to Fe 1 abundances 
and see this as validation that LTE assumption s still hold when 
studyi ng this type of star. In contrast, work bv lShchukina et all 
d2005l) find higher correction values of ~0.9 dex and ~0.6 dex, 
depending on whether 3D or ID atmospheres are used. The dif- 
ference in their conclusions is driven principally by the differ- 
ent relative import ance of collisional an d radiative transitions in 
their calculations. IGratton et alJ (JT999) have relatively stronger 
collisional transitions, and a s a r esult find smaller deviations 
from LTE. Shchu kina et alJ (120051) include no collisions with 
neutral hydrogen. We return to this important point below, but 
for now it illustrates that much work still needs to be done in 
this field before we can be certain of the impact of NLTE. 

2.2. The coupling of the radiation field and level populations 

Many factors have to be taken into account when computing ra- 
diative transfer in NLTE. This leads to a complicated situation 
where, for example, we have to solve population equations and 
radiative transfer equations simultaneously. This is due to the 
level populations and the radiation field being coupled, a fact ig- 
nored in LTE calculations. There are large uncertainties in NLTE 
calculations because of the lack of complete information on the 
rates of collisional and radiative transitions between energy lev- 
els for a given element in all its important ionization states. This 
is especially true for larger atoms which have a greater number 
of energy levels, as is the case for Fe. 

To solve NLTE problems, a system of rate equations is 
needed that describes fully the populations of each level within 
the atom under study. Statistical equilibrium is invoked, i.e. the 
radiation fields and the level populations are constant w ith time. 
The f ormulation of the problem is well described in Mihalas 
dl978l) . from which the following equations are taken. The pop- 



A. Hosford et al.: Lithium abundances of halo dwarfs based on excitation temperatures. II. NLTE 



3 



ulation of level i is the sum of all the processes that populate the 
level minus the processes that depopulate it, such that: 



dw,- 
df 



1*> 



N 



(1) 



where n, and n ; are the populations of the levels i and j respec- 
tively, N is the total number of levels, including continua, and 
Pji and Pij are the rates of transitions into and out of the level i. 
The rates are given by: 



Pij — Ajj + BijJ V(j + Ci 



(2) 



where Ay, By, and Cy are the Einstein coefficients for sponta- 
neous, radiative and collisional excitation respectively, vo is the 
frequency of the transition between levels i and j, and J Vo is the 
mean intensity averaged over the line profile. It is usually the 
case that the radiative rates dominate over the collisional ones 
at optical depths 75000 < 1 implying that LTE assumptions no 
longer hold in general. The rate equations depend on the mean 
intensities over the relevant frequencies, J V(I , meaning that the 
level populations depend on the radiation field. Conversely the 
radiation field depends on the level populations through the ra- 
diative transfer equation. This is seen by examining the simple 
problem of the two-level atom where the radiative transfer equa- 
tion is given by: 



(3) 



where 7 y is the intensity at frequency v, /j. is the cosine of the 
viewing angle, t y is the optical depth and S Y is the line source 
function, such that: 



Sl = 



<2hv 2 



n i8j 
njgi 



- 1 



(4) 



Here h is Planck's constant, c is the speed of light, g t and gj are 
the statistical weights of the levels i and j respectively, and tii/rij 
gives the ratio of the populations of the levels i and j calculated 
using the rate equation, Eq. (1). This form of the source function 
strictly speaking holds under the assumption of complete fre- 
quency redistribution but still illustrates the problem of having 
to solve the two sets of equations simultaneously. 

2.3. Transition rates 

For the calculation of the level populations, through Eq. 1, radia- 
tive and collisional rates are required. 

For the radiative rates, the bound-bound transition proba- 
bilities and photoionization cross sections are needed for all 
levels of the atom in all significant ionization states. Two of 
the larg er projects pro viding values for these are the Opacit y 
Project dSeatonl[T987h and the IRON project (lBautistalll99^) . 
For Fe, the Opacity Project fin ds typically a > 10 % uncertainty 
for their photoionization data (ISeaton et al.lll994h . The Bautista 
photoionization values, which are larger th an those previ ously 
used, lead to increased photoionization rates ( Asplund 2005]) and 
hence to lower abundances as overionization becomes more ef- 
ficient. 

For the collisional data, large uncertainties still exist. The 
two main types of collisions that affect the line profile are those 
with electrons and neutral hydrogen. Coupling of all levels in the 
Fe model atom occurs due to these types of collisions, especially 
in the atmospheres of cool stars where electrons and neutral H 



are believed to be the dominan t perturbers. A simple calcula- 
tion, like that in lAsplundl(l2005h . shows that H I collisions dom- 
inate over electron collisions in thermalizing processes in metal- 
poor stars and are therefore important in calculations of line 
profiles. For co llisions with neutral h ydrogen, the approximate 
fo rmulation of | Drawinl (119681 1 1969b is used as implemented 
bv lSteenbock & Holwegerl(ll984l) . However, through laboratory 
testing and quantum calculations of collisions with atoms such 
as Li and Na, it has been shown that Drawin's formula does not 
produce the correct order of magnitude result for H I collisional 
cross-sections. In some cases, where comparisons with experi- 
mental data or theoretical results can be made, the Drawin recipe 
overestim ates the cross-sections by one to six o rders of magni- 
tude (e.g. lFleck et alJll991tlBlirklem et al.ll2003l). Correction s to 
the Drawin cross-sections are suggested by Lambert ( 1 1993b to 
compensate for these differences. 

Due to the uncertainties in the magnitude of the H collisions, 
the Drawin cross-sections are scaled with a factor Sh- There are 
differ ent schools of th ought on how to deal with this parame- 
ter. ICollet etaf] (120051) treat it as a free parameter in their work, 
adopting values of Sh = 0.001 and 1 and test the effect this has on 
their results. Higher values of Sh correspond to more collisions 
and hence more LTE-like conditions. Their main aim, however, 
was to test not the efficiency of H c ollisions but the e ffects of 
line-blocking on the NLTE problem. iKorn et alJ J200l make it 
one of their aims to constrain Sh- To do this, they ensure ioniza- 
tion equilibrium between Fe 1 and Fe 11 using the log g derived 
from hipparcos parallax and T e ff from H lines. In doing this, they 
find that a value of Sh = 3 holds for a group of local metal- 
poor stars. This apparently contradicts the statem ent above that 
Draw in's formula overestimates the cross-sections. Gratto n et al.l 
(1999) use Sh = 30. This value was constrained by increasing Sh 
until spectral features of several elements, i.e. Fe, O, Na and Mg, 
of RR Lyrae stars all gave the same abund ance. With such ele- 
vated collisional rates, Gratto net al.l l[l999) not surprisingly find 
results very close to LTE, i.e. they find very small NLTE correc- 
tions. 

Collisions with neutral hydrogen and electrons are impor- 
tant not only in coupling bound states to each other, but also 
in coupling the whole system to the continuum i.e. to the Fen 
ground state (and potentially excited states). This is especially 
true when considering the high excitation levels. These levels 
are more readily collisionally ionised than lower levels, and are 
also coupled to each other by low energy (infrared) transitions, 
therefore thermalization of the levels occurs which drives the 
populations more towards LTE values. It is therefore important 
to have a model atom th at includes as many of the higher terms 
of the atom as possible dKornl2008h . although it is not necessary 
to include all individual levels. We return to this point in Sect. 4. 
We describe the model atom and calculations next before mov- 
ing on to the results. 

2.4. The model Fe atom 



The F e model adopted for this work is that of ICollet etafl 
(2005), whic h is an up dated version of the model atom of 
iThevenin & Idi art ( 1999). The atom includes 334 levels of Fei 
with the highest level at 6.91 eV. For comparison, the first ioniza- 
tion energy is 7.78 eV and the NIST database lists 493 Fe 1 levels. 
Many of the highest levels are not included in our model; due in 
part to computational limitations i.e. the more complicated the 
model, the greater the computer power and time needed to com- 
plete the computations, and because of lack of important infor- 
mation, e.g. photoionization cross sections. We report below on 
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the effects the missing upper levels have on the corrections and 
try to quantify their importance in the NLTE calculations. The 
model also includes 189 levels of Fen with the highest level at 
16.5 eV, and the ground level of Fein. For comparison, the sec- 
ond ionization energy is 16.5 eV, and the NIST database lists 578 
Fen levels. This model configuration leads to the possibility of 
3466 bound-bound radiative transitions in the Fe i system, 3440 
in the Fen system, and 523 bound-free transitions. We run the 
calculations with the whole model, but present results only for 
the lines that are measured in our program stars. 

Oscilla t or str engt hs for the Fei line s are taken from 
Nave et al] (1 1994b a ndlKurucz & Belli dl995l). w hilst values from 
Fuhr etal] d 19881) . iHirata & Horaguchil d 19951) . and iThevenrril 
dl989l [19901) were used for the Fen lines. The photoionization 
cross-sections are taken from the IRON Project (Bautista 1997). 
ICollet et alj j2005) smoothed these cross-sections so as to mini- 
mize the number of wavelength points to speed up the computa- 
tional processes. 

Collisional excitations by electrons are incorpora ted through 
the van Regemorter formula (van Reg emortenl 19621) and cross- 
sections for collisi onal i onization by electrons are calculated 
by the methods of ICoxl d2000j). In the case of H collisions, 
the appro ximate description of iDrawinl d 19681 [l969), as imple- 
mented bylSteenbock & Holw eger (1984) with the correction of 
Lambert (1993) and multiplied by Sh, has been used. As we 
do not intend to constrain Sh, we treat it as a free parameter 
and adopt values of (no neutral H collisions), 0.001 and 1 
(Drawin's prescription). This allows us to assess the importance 
of H collisions on the NLTE corrections. For all calculations, 
the oscillator strength value, fij, has been set to a minimum of 
10~ 3 when there is no reliable data or the / value for a given line 
is below this minimum. This minimum is set as the scaling be- 
tween the cross-sect ions and the / value breaks down for weak 
and forbidden lines dLamberjll993h . 



2.5. The model atmospheres 

In this work, we have adopted plane-parrallel MARCS models. 
These model s are used, rat h er tha n the Kurucz 1996 models as 
was done in iHosford et al. (2009), as MULTI needs a specific 
format for its input, this is provided by the M ARCS, details of 
which can be found in lAsplund et al.l dl997l) . 3D models lead 
to an even steeper temperature grad ient, and hence cooler tem- 
peratures in the line forming region (Asplund 2005), but the use 
of these more sophisticated models is beyond the scope of this 
work. 



2.6. Radiative transfer code 

The NLTE code used to produce Fe line profil es and equivalen t 
widths (W^) is a modified version of MULTI dCarlssonl fl986). 
This is a multi-level radiative transfer program for solving the 
statistical equilibrium and radiative transfer equations. The code 
we adopted is a versio n modified by R. C ollet to include the ef- 
fects of line-blocking dCollet et al]|2005l) . To do this, they sam- 
pled metal line opacities for 9000 wavelength points between 
1000 A and 20000 A and added them to the standard background 
continuous opacities. They found that, for metal-poor stars, the 
difference between NLTE Fe abundances derived from Fe i lines 
excluding and including line-blocking by metals in the NLTE 
calculations is of the order of 0.02 dex or less. 



3. NLTE calculations 

For this work, we ha ve chosen six of our original program stars 
dHosford et al]|200"9l) that approximately represent the limits of 
our physical parameters, i.e. one of the more metal-rich, one of 
the less metal-rich, one of the hotter, one of the cooler etc. Table 
[TJ indicates the stellar parameters for which model atmospheres 
were created. For the HD stars hipparcos gravities were used. For 
the other three stars, lower an d upper Umits on log g are given 
by theoretical isochrones (see IHosford et al.ll2009l) . In the case 
of LP8 15-43, there is uncertainty as to whether it is just above 
or just below the main-sequence turnoff. The final temperatures 
are interpolated between these values using a final log g that 
represents the star at 12.5 Gyr (Table|2]). This study is primarily 
concerned with the formation of Fe i lines. 



Table 1. Physical parameters for the atmospheric models used 
in this work 



Star 




[Fe/H] 


l°g<? 






(K) 


(dex) 


(dex) 


(Kms- 1 ) 


HD 140283 


5769 


-2.54 


3.73 


1.5 


HD84937 


6168 


-2.34 


3.98 


1.3 


HD74000 


6070 


-2.20 


4.03 


1.2 


BD+26° 2621 


6225 


-2.68 


4.47 


1.2 


BD+26° 2621 


6241 


-2.67 


4.51 


1.2 


CD-33° 1173 


6380 


-2.94 


4.41 


1.5 


CD-33° 1173 


6391 


-2.94 


4.47 


1.5 


LP815-43 (SGB) 


6383 


-2.71 


3.80 


1.4 


LP815-43 (SGB) 


6409 


-2.68 


3.91 


1.4 


LP815-43 (MS) 


6515 


-2.62 


4.35 


1.4 


LP815-43 (MS) 


6534 


-2.61 


4.42 


1.4 



In Fig.[TJ we present the departure coefficients, bj = nj/nf^, 
for the lower (left hand side) and upper (right hand side) levels 
of all lines we have measured in the star HD 140283 in paper 
I, calculated for three Sh values. The two sets of lines in each 
plot, coloured, red and blue, represent levels that fall above and 
below the midpoint of our excitation energy range, i.e. 1.83 eV 
where our highest lower level of the transition is at 3.65 eV, and 
5.61 eV where our highest upper transition level is at 6.87 eV. 
This is done to better visualise the effects of NLTE on different 
levels of the atom. We see that in all cases the Fei levels are 
under-populated compared to LTE at T5000 <1- This is primar- 
ily due to the effects of overionization where J v > B v for lines 
formed from the levels of the atom at around ^ ~ 4 eV below 
the continuum, due to the UV photons having energies « 3-4 
eV. This causes all levels of the atom to become greatly depop- 
ulated, as can be seen from the blue lines. The coupling of the 
higher levels through collisions and of the lower levels through 
the large number of strong lines sharing upper levels implies that 
relative to one another the Fe 1 level populations approximately 
follow the Boltzmann distribution. Because of photoionization, 
the Saha equilibrium between Fe 1 and Fe n is not fulfilled how- 
ever and the departure coefficients of Fei levels are less than 
unity. In deeper levels of the atmosphere, this leads to both up- 
per and lower levels of a transition being equally affected by 
the above phenomena (Fig[TJ- right hand side). For this reason, 
the source functions for lines forming at these depths are rel- 
atively unaffected in this region, as S * {b uvva lb\ ovll£S )B Y , and 
follow a Planckian form (Fig. [2]— right hand panel). The com- 
bined effect of the above processes, i.e. depopulation and rel- 
atively unaffected source functions, leads to a smaller Wa and 
thus weaker lines, and increased abundances compared to the 



A. Hosford et al.: Lithium abundances of halo dwarfs based on excitation temperatures. II. NLTE 



LTE case. For stronger lines, forming further out in the atmo- 
sphere, there is a divergence between b upper and bi ower and the 
source function thus diverges from the Planck function (Fig. [2]— 
left hand panel). In the case where S' v < B v , the source function 
compensates slightly for the loss of opacity leading to smaller 
NLTE corrections, the opposite being true for S' v > B v . We see 
that for the lower level of the weaker line considered in the fig- 
ure has J v > By, whilst S [ ~ B v , which leads to overionization 
of that level and greater departures than the stronger line and 
greater NLTE abundance corrections. 

The effect of H collisions is in general to reduce the spread 
of departure coefficients and drive populations towards LTE val- 
ues. This reduction in the spread of departure coefficients comes 
from the coupling of bound states. The increase of H collisions 
gradually reduces the departures from LTE through the atmo- 
sphere as shown in Fig.[TJ with an increasing Sh the slope in the 
departure coefficient profile becomes shallower. In Fig.[TJit is in- 
teresting to see that the rise in b t at around 7-5000 ~ -2.5 for the 
levels below 1.83 eV becomes smaller with increasing Sh- This 
could in fact mean an increase in NLTE departures for some lev- 
els for increasing S h, rather than H collidions driving conditions 
towards LTE which is normally the case. This rise is most likely 
caused by increased recombination in the upper (infrared) levels 
followed by a cascade of electrons down to lower levels. Exactly 
how this is affected by the increase in Sh is not yet known and 
requires further study. 

The decrease in level population at 75000 < 1 causes a drop 
in opacity for all lines. As a result of this, the lines form deeper 
in the atmosphere than in LTE. In Fig. [3] we clearly see this 
effect, where we show the continuum optical depth T5000 at 
which the line optical depth r v = 2/3. We also see that there 
is an increasingly large logarithmic optical depth difference, 
A log T5ooo(i"v = 2/3), between the formation of weak lines in 
NLTE and LTE, up to m 50 mA, after which the difference be- 
comes constant. With a decrease in opacity compared to LTE, 
there needs to be an increase of abundance to match the equiv- 
alent width of a given line in NLTE. Opacity is not the only 
variable affected by NLTE, the source function can also be af- 
fected. However, it is the dominant force in driving the NLTE 
departures within the Fe atom. In Fig.|4j we plot the abundance 
correction versus equivalent width for the star HD140283. We 
see that there is a positive correction for the different values of 
Sh- There is a clear trend with equivalent width. It is how this 
translates to trends with excitation energy % that will affect Teg: 
if the abundance corrections only shifted the mean abundance 
without depending on^ then the derived r e tf would not change. 

Through Fig.[T]to Fig.|4]the general effects of NLTE on line 
formation can be seen. The depletion of level populations (Fig. 
Q]) leads to a lower opacity and shifts the depth of formation to 
deeper levels (Fig. [3]). This also means that a higher abundance 
is needed within NLTE, leading to positive abundance correc- 
tions (Fig. |4|. However, there is a competing effect in some 
cases where the source function deviates from the Planck func- 
tion (Fig. [2j, which, in the case of the strong lines, compensates 
for the level depletion and decreases the abundance correction, 
as is seen in Fig. [4] 

4. NLTE abundance corrections - deriving, testing, 
applying 

In order to determine a new for a star, we first need to calcu- 
late NLTE corrections for the LTE abundances derived in Paper 
I. Abundance corrections of the form Anlte.multi - ^lte.multi 
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Fig. 3. The depth of formation of Fe 1 lines, with no H collisions, 
on the log T5000 scale for lines of different equivalent width, in 
LTE (filled circles) and NLTE (open circles), for HD140283. 
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Fig. 4. Abundance correction versus equivalent width for the 
lines measured in the star HD 140283 for Sh = 1 (filled trian- 
gles), 0.001 (filled circles) and (open circles). 

are calculated and applied to the LTE abundances from Paper I 
to generate NLTE abundances on the same scale as that paper, 
rather than using solely the new NLTE analysis. This procedure 
is used so as to tie this work to the previous results, thus allowing 
the limitations of the LTE assumptions in that work to be seen. 
To do this, a grid of MULTI results for a range of abundances 
is created with increments of 0.02 dex. The abundance values 
covered by this grid depend on the spread of abundances from 
individual lines in each star. MULTI gives an LTE and NLTE 
equivalent width for each abundance in this grid. A first step is 
to determine what Wlte from the M ULTI grid correspon ds to the 
LTE abundance derived in Paper I dHosford et alj|2009l) . This is 
done for all Fe lines that are measured in the star. The NLTE 
abundance inferred for a line is the abundance that corresponds 
to this W,\ within the grid of MULTI NLTE results. The correc- 
tion is then calculated as AA(Fe) = A(Fc)nlte - A(Fc)lte- Fig. [5] 
shows the corrections for the star HD 140283 calculated for the 
three different Sh values: Sh = 0, 0.001 and 1. We see a trend 
in the abundance correction with^, where we have values, from 
least square fits, of: 

A(Fe) = 0.490(+0.012) + 0.0216(+0.0055)^; for S H = (5) 
A(Fe) = 0.490(+0.012) + 0.0207(+0.0054)^; for S H = 0.001(6) 
A(Fe) = 0.244(+0.072) + 0.0178(+0.0032)^; for S H = 1 (7) 
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The non-zero coefficient of x implies that a T e g correction is 
needed. The values for Sh = and Sh = 0.001 are very similar 
and imply that corrections for these two values will be very 
similar. We therefore decided that corrections for only Sh = 
and 1 would be calculated, Sh = representing the maximal 
NLTE corrections and Sh = 1 representing the full Drawinian 
magnitude of neutral H collisions. 
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Fig. 5. Abundance correction versus % for Sh = 1 (filled trian- 
gles), 0.001 (filled circles) and (open circles) for the measured 
lines in the star HD140283. Least-squares fits are shown to guide 
the eye. 

To test the corrections, we compared synthetic profiles from 
the NLTE abundance with the observed profile, and compared 
measured W/s with NLTE W/s from MULTI, obtained from 
an abundance given by Alte + AA. The synthetic profiles are 
convolved with a Gaussian whose width is allowed to vary from 
line to line. This represents the macroturbulent and instrumen- 
tal broadening, the latter calculated by fitting Gaussian profiles 
to ThAr lines in IRAF and found to be ~ 100 mA. We found 
that the profiles match the observed line reasonably well, and 
that measured and MULTI calculated W/s are comparable, with 
a standard deviation of 2.3 mA. This gives us confidence that 
the corrections are realistic within the framework of the atomic 
model used. These corrections were then applied to the WIDTH6 
LTE abundances used in Paper I and new plots of ^ versus A(Fe) 
were plotted. We then nulled trends in this plot to constrain 
r e ff(NLTE) by recalculating the LTE abundances using the ra- 
diative tran sfer program WIDTH 6 (Kurucz & Furenlid 1978) 
exactly as in lHosford et al.l (120091) and reapplying the NLTE cor- 
rections, derived here from MULTI for the original LTE param- 
eters. 

As noted in Sect. 12. 31 it can be important to include the high- 
est levels of the atom in the calculations. It is not necessary to 
include each individual level however, and it is possible t o use 
super levels that represent groups of closely spaced levels dKornl 
2008). To test the effect of these upper levels, we took the ap- 
proach of giving the top 0.5 eV of levels in our atomic model an 
Sh = 2 whilst the rest of the levels had Sh = 1. We have done 
this for three situations; A) increasing Sh for just the bound- 
bound transitions rates, B) increasing Sh for just the bound-free 
rates and C) increasing Sh for both the bound-bound and bound- 
free. We discuss here only the case of the bound-free rates as 
it is only these rates that have an effect, edging the populations 
towards LTE values. Changing the bound-free rates not only af- 
fects the higher levels but translates through all lower ones. In 
fact it is the lower half of the atomic model that is affected by a 



greater amount; further investigation into reasons for this effect 
are discussed in Sect. 16. II The result can be seen in Fig.|6]where 
we plot a level with x — 0-96 eV and one of the higher levels, x 
= 3.30 eV, from our atomic model. Fig.|7]shows the abundance 
correction against^ for the increased Sh value of the upper lev- 
els and for a pure Sh = 1 situation. Comparing the differences 
in abundance correction between Sh = 1, and Sh = 1 with Sh 
= 2 on the upper levels we see a mean difference (AA(Fe)s H =i+2 
- AA(Fe)s H= i) of -0.031 dex for x = - 2 eV and -0.028 dex 
for^ > 2 eV, for the star HD140283. These effects equate to a 
5 K increase in T e s compared to Sh = 1. It is then clear that the 
upper levels have a slight effect on the final temperatures, and 
induce a slightly larger NLTE correction. However, in the case 
of this study, where random errors are of order ~ 80 K, they will 
not make a significant effect. 
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Fig. 6. Dashed line: The effects of increasing the Sh value to 
2 for the top 0.5 eV of levels in the atomic model, for bound- 
free transitions from a low level (top panel) and a higher level 
(bottom panel). The remaining levels have Sh = 1. Solid line: 
Sh = 1 for all levels. 



5. Results 

Fe abundance corrections for the stars in Table Q~]have been cal- 
culated and new temperatures have been derived using the exci- 
tation energy technique, as in Paper I but with the NLTE correc- 
tions applied as described in Sect. [4] Table |2]lists the new NLTE 
r e ff's and Ar e fr, such that Ar e fr = T e ff(NLTE) - T e g(LTE), for 
the selection of stars. 

For this work, all the other par ameters, viz. logj?, [F e/H] and 
were kept at the values found in lHosford et alJ c2009). Our aim 
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Table 2. Final T e ff and A(Li) for the selection of stars in this study 
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Fig. 7. Comparison between the abundance correction versus ex- 
citation energy for the star HD 140283 using Sh = 2 for the top 
0.5 eV of levels of Fe i in the atomic model, and Sh = 1 for the 
remaining levels (crosses) and pure Sh = 1 (open triangles). 

here, as it was in Hosford et al. (2009), is to narrow down the 
zero point of the temperature scale by quantifying the system- 
atic errors, albeit at the expense of having larger star to star ran- 
dom errors. Contributions to the errors come from adopted grav- 
ity, the nulling procedure in determining the T e ff, and smaller 
contributions from the error in microturbulence, errors in the 
age, metallicity and initial temperature, r p h t, when determin- 
ing isochronal gravities. In relation to the gravities, the three 
HD stars had gravities derived using hipparcos parallaxes, and 
their errors are a reflection of errors propagating through this 
calculation, whilst for the remaining stars isochrones were used. 
The isochronal gravities are sensitive to age, with a 1 Gyr differ- 
ence leading to a change of ~ 0.03 dex for main sequence (MS) 
stars and ~ 0.06 dex for sub-giant (SGB) stars. This equates 
to a change in T x of 12 K and 24 K respectively. These errors 
are based on LTE sensitivities, as are other errors quoted below. 
There is also a dependence on the initial temperature, a photo- 
metric temperature from Ryan et al (1999), used to determine 
the isochronal gravity. A +100 K difference leads to +0.06 dex 
and -0.06 dex for MS and SGB stars respectively. This equates 
to ± 24 K in T x which shows, importantly, that T x is only weakly 
dependent on the initial photometric temperature. Contributions 
to T x is also sensitive to microturbulence, for which an error of 
~ 0. 1 km s _1 equates to an error of » 60 K. 

In the nulling procedure any trends between [Fe/H] and x 
are removed. Due to the range in line to line Fe abundances for 
a particular star, there is a statistical error in the trend which is 
of order <x = 0.01 1 dex per eV, which equates to « 40 K - 100 
K depending on the star under study. This error also contains 
the random line-to-line errors due to equivalent width, gf, and 



damping values. The final error in Table|2]is then a conflation 
of this statistical error and the errors from Aage = 1 Gyr, A£ = 
0.1 kms" 1 , A [Fe/H] = 0.05 and A7p hot = 100 K. 

Th ese new T e g values and equivalent widths from lRvan et al.l 
(1999) were then used to calculate new Li by interpolating 
within a grid of equivalent wi dth versus abundan ce for different 
T e g. This grid was taken from lRvan et al.l (fl996a). 

6. Discussion 

6.1. The r eff scale 

With the addition of the NLTE corrections, we see in Table [2] 
that the re is, for the most par t, an increase in r e jj from the LTE 

's of lHosford et all {2009), for both cases of Sh- The only ex- 
ception is CD-33 1 173 in the Sh = case, for which there is a 
93 K decrease. We return to this star below. The r e jf corrections 
we have derived average 59 K for Sh = and 73 K for Sh = 
1 (treating LP8 15-43 as one datum, not two). For Sh = the 
temperature corrections tend to increase at cooler temperatures, 
whilst the tendency is weaker or opposite for Sh =1, i-e- correc- 
tions increase at higher temperatures (obviously the gravity and 
metallicity of the stars also affects their NLTE corrections, but 
nevertheless we find it intsructive to consider temperature as one 
useful discriminating variable). This gives rise to a change in 
the difference Ar e ff,s H =o _ ATeff,s H =i w i m temperature, with this 
quantity being negative for the two hottest stars, CD-33°1173 
and LP815-43. The switch over from Sh = having the larger 
correction to Sh = 1 having the larger correction is at around 
r e ff « 6200 K . Further testing has shown that this is not a 
random error and is clearly something to investigate further in 
the future. This is further shown by Fig. [8] where the abundance 
correction versus^ for the stars CD-33°1173 and LP815-43 
(SGB) are plotted. It is seen that for LP815-43 (SGB), increas- 
ing Sh has a larger effect on the lower excitation lines than for 
higher ones. This has induced a trend of abundance with^- larger 
than that of the Sh = case. This in turn leads to a larger tem- 
perature correction for Sh = 1 than for Sh = 0. The reason for 
this effect is still uncertain. 

To investigate this behaviour further, the test of increasing 
the Sh value of the upper levels, as done on HD 140283 in Sect 4., 
has also been performed on LP815-43 for the MS and SGB pa- 
rameters. This has shown that the effect of collisions with neutral 
H are indeed larger for the lower levels of the atom, and that this 
effect is larger for LP8 15-43 (MS), which is the hottest star. This 
indicates that there is a temperature dependence, i.e. the differ- 
ence between the mean difference (AA(Fe)s H =i+2 - AA(Fe)s H= i) 
(where Sh = 1 + 2 indicates the scenario of having Sh = 2 for 
the top 0.5 eV worth of levels) for the levels with y < 2 eV 
and those with^ > 2 eV is greater for the hotter star, LP815-43 
(MS). However, when performing this test on LP8 15-43 (SGB), 
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which has a similar log g to HD 140283 whilst still being hotter, 
the effect is not as great as for HD140283. This shows that there 
is some gravity dependence on the neutral H collisions along 
with the temperature dependence i.e. the gravity indirectly af- 
fects the collisional rates, by impacting on the number density 
of hydrogen atoms at a given optical depth. Fig. [8] along with 
Fig. [5] clearly show that NLTE has varying star to star effects, 
i.e. from the similar effects at different Sh values in HD 140283 
(Fig.|5j, to the differing effects in CD-33°1173 and LP815-43 
(SGB) (Fig. [8). The range of AT e g values, and the negative value 
for CD-33 1 173, shows the intricacies of the NLTE process, and 
that generalisations are not easily made when identifying the ef- 
fects of NLTE on temperatures determined by the excitation en- 
ergy method. For the purposes of this paper, which is concerned 
with the effective temperatures in the context of the available 
NLTE model, it is appropriate to acknowledge these NLTE ef- 
fects and to move ahead to use them in the study of the Li prob- 
lem, whilst still recognising that much work remains before we 
approach a complete description of the Fe atom. 

Although we discussed the possibility that the extreme (neg- 
ative) Ar e ff correction for CD-33°1173 is due to corrrections 
being temperature-dependent, this unusual case may be in part 
due to the fact that only a subset of the original lines mea- 
sured is available through the NLTE atomic model. The atomic 
model does not contain every level of the Fe atom and there- 
fore some transitions are not present in the calculations. This 
means that not every line measured for a given star is present 
in the calculations and leads to a trend being introduced in the 
^-abundance plot prior to the trend induced by the NLTE cor- 
rections. This is because the original nulling of the ^-abundance 
plot was achieved with a greater number of points. CD-33° 1 173 
has the least lines available from the atomic model used with 
MULTI, however, there is no distinct trend between AT e f{ and 
the number of lines available for each star, and after testing we 
found that the effect of the subset, i.e. the measured lines that are 
available with our atomic model, is to increase the LTE temper- 
ature. This implies that the decrease in T e g for this star is most 
likely due to NLTE effects. Although there is no obvious corre- 
lation between the number of lines available and the tempera- 
ture correction, this emphasises the need for a complete atomic 
model. This is especially true when considering the abundance 
of individual lines, as in the excitation technique used in this 
work. 

As i n Paper I, we ha v e compared our T e g values with 
those of iRvan et all (Il999h . IMelendez & Ramfrezl d2004 . and 
lAsplund et al.l ( 2006 ). Fig. [9] presents these c omparisons- 
Comp aring against the photometric temperatures of Rya n et aTl 
( 1999) for five stars in common, we see that our new T e g scale is 
hotter by an average of 132 K, with a minimum and maximum 
of 43 K and 21 1 K respectively for an S H = 0. Recall that S H = 
corresponds to the maximal NLTE effect, i.e. no collisions with 
the hydrogen, for the model atom we have adopted. For Sh = 1, 
our scale is hotter by an average of 162 K, with a minimum and 
maximum of 101 K and 267 K respectively. 

We have three stars in common with Melen dez & Ramfre A 
(2004). Their temperatures are hotter than the ones we derived 
here by 196 K on average for Sh = with a minimum and 
maximum difference of 27 K and 381 K respectively, and by 
193 K on average for Sh = 1, with a minimum and maxi- 
mum difference of 84 K and 247 K respectively. Therefore, 
even with NLTE corrections we still c annot achieve the high 
T e g of the IMelendez & Ramfrezl (120041) study. It has however 
been noted (Melende z 2009 - private communication) that the 
IMelendez & Ramfrezl ([2004) temperatures suffer from system- 
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Fig. 8. Abudance correction versus x f° r (top panel) 
CD-33°1173 and (bottom panel) LP815-43 for S H = 
(filled circles) and Sh = 1 (open circles). 



atic errors due a imperfect calibration of the bolometric correc- 
tion for the choice of photometric bands used. This led to an 
inaccurate zero point and hotter Tiff's than most other studies. 
The revision of their temperature scale is not yet available and 
comparisons to their new T^'s is not possible at this time. 

Fin ally, we have three stars in common with Asp lund et al.l 
(120061) . Using Sh = we obtain temperatu res for two of the stars 
that are hotter than lAsplund et al.l d2006h by 97 K and 151 K. 
The third star is CD-33°1173, for which we calculated a neg- 
ative temperature correction, and which is cooler in our study 
by 97 K. The temperatures f or all three stars are hotter in our 
study than in Aspl und et al.l (120061) when using Sh = 1. Here 
the average di fference is 110 K, va lues ranging from 37 K to 
207 K. If the lAsplund et al. Uiaiemperatures are affected 
by NLTE, as stated bv iBarkleml d2007h who expects a 100 K 
increase in Balmer line temperatures, this would bring the T e g 
scales back into agreement. Another problem facing the Balmer 
line method is the effects of granu lation, due to convection, on 
the line wings dLudwig et al.l2 009). It has been found (Bonifacio 
- private communication) that inclusion of these effects would 
increase the effective temperatures derived with this method. In 
particular a value of r e ff = 6578 K has been found for the star 
LP815-43. Although this is 176 K hotter than our result for the 
SGB case with Sh = 0, i.e. T e g = 6402 K, it is in good agreement 
with the values T e ff = 6522 K (Sh = 1) for the SGB case and T e g 
= 6551 K (S H = 0) or T e g = 6607 K (S H = 1) for the MS case, 
calculated in this work. 
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Fig. 9. T e ff comparis on between this work andlRyan et al.l (j 1999. 
R99) (filled ci rcles). IMelendez & Ramirez! (|2004] MR04) (filled 
triangles) and lAsplund et alJ (120061 A06) (filled squares), for 
(top panel) Sh = and (bottom panel) Sh = 1. Dashed lines 
are + 100 K limits. 

6.2. Lithium abundances 

We now address the new Li abundances and their effect on the 
lithium problem. We see that the introduction of NLTE correc- 
tions to the r e ff scale has led to temperatures that are of order 
100 K hotter t han LTE temperature s cales, with the obvious ex- 
ception of the Mele ndez & Ramfrez] (|2004) scale. This will then 
lead to an increase in the mean lithium abundance. Table[2]lists 
A(Li) for the new temperatures. With these new r e ff's, we calcu- 
late a mean Li abundance of A(Li) = 2.19 dex with a scatter of 
0.072 dex when using Sh = 0, and A(Li) = 2.21 dex with a scatter 
of 0.058 dex for the Sh = 1 case. Consistent with the temperature 
increase, these v alues are higher th an those found by other stud- 
ies, in particular lSpite et aU (ll99o) . who found a value of A(Li) 
= 2.08 (±0.08) dex using a similar iron excitation energy tech- 
nique but without the NLTE corrections, Bonifa cio et alJ d2007l) 
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Fig. 10. Lithium abundance versus [Fe/H] (top panel) and T^g 
(bottom panel) for Sh = (filled circles) and 1 (open circles). 

with A(Li) = 2.10 (+0.09) using a Balmer line wing temperature 
scale, and A(Li) = 2. 16 dex or A(Li) — 2.10 depending on the 
evolutionary state from iHosford et alJ (120091) . The NLTE cor- 
rections have moved the mean Li abundance closer to, but not 
consistent with, the WMAP value of A(Li) = 2.72 dex, and thus 
stil l leaves the lithium problem unsolved. It is noted that even 
the IMelendez & Ramfrezl (120041) scale, whilst bringing the ob- 
served and theoretical Li abundances closer, still failed to solve 
the lithium problem. 

Fig.[10]shows the lithium abundances versus [Fe/H] and 7T e ff, 
least squares fits have been performed for both sets of data. In the 
fit to metallicity we get the values: 

A (Li) = 2.55(±0.31) + 0.137(±0.119)[Fe/H] (8) 
for Sh = and 

A(Li) = 2.04(+0.26) - 0.067(+0.102)[Fe/H] (9) 

for Sh = 1. For Sh = 0, we see a sim ilar values to the coeffi- 
cient of metallicity as iRvan et alJ (Il999l) . whilst for Sh = 1, we 
have a value that is about half the size, and has a negative slope. 
However, our errors are much larger, due to the large random 
errors and small sample of stars, and therefore no statistically 
relevant trend can be deduced. For 7T e ff, we get the equations: 

A(Li) = 2.77(+0.91) - 0.00009(+0.00015)T eff (10) 

for Sh = and 

A (Li) = 2.09(+0.63) + 0.00002(+0.00010)T eff (11) 

for Sh = 1 ■ Here we see no statistically relevant trend with 
for either Sh value. 
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We also perform the fit as described by iRvan et aD ([2000), 
such that: 

Li/H - a + b Fe/FeQ (12) 

where a measures the primordial abundance of Li and b' is a 
probe of galactic chemical evolution. For this fit, we obtain the 
primordial values of 7 Li/H = (1.47 ± 0.27)xl0" 10 forS H = and 
7 Li/H = (1.80 + 0.16)xlQ- 1Q for S H = 1. Both of th ese values are 
far fro m the high value of 7 Li/H = 5.24+° 10~ 10 (ICvburt et al.1 
2008)from WMAP and BBN. 

We see then that the addition of NLTE corrections has led to 
an increase in T e ff for most stars. This equates to an increase in 
A(Li) but it is still not high enough to reconcile th e lithium prob- 
lem. Through the efforts of Hosf ordet al.l (|2009) and this study 
it is safe to conclude that systematic errors in the metal-poor r e ff 
scale are almost certainly not large enough to be the source of 
the A(Li) discrepancy between observation and WMAP + BBN 
predictions. This outcome lends strength to other possible expla- 
nations, such as processing in the stars, e.g. diffusion, processing 
in earlier generations of stars, and/or different BBN networks, 
or more exotic solutions requiring physics beyond the standard 
model. 

It should be noted that while we have computed Fe lines in 
NLTE to constrain the temperature, our Li abundances are cal- 
culated from a grid of abundance versus equivalent width that 
was constructed under the assumptions of LTE, see Paper I for 
details. Several studies of the effects of NLTE Li line forma- 
tion have been cond ucte d. Two of t hese s tudies are those of 
ICarlsson et all (11994 and lLind et aD d2009t) : they find Li abun- 
dance corrections of ~ +0.013 — f-0.020 dex and ~ +0.01 - 
+0.03 dex respectively for the temperature, log g and [Fe/H] 
range in this study. Due to the very small size of these correc- 
tions we find the use of the LTE grid, combined with our NLTE 
effective temperatures, to be acceptable in determining Li abun- 
dances, and that the introduction of NLTE Li abundances will 
not significantly aid in solving the lithium problem. 

6.3. Implications of NLTE calculations for ionization balance 
and S H 

Having discussed the effects of NLTE corrections on the T e g 
scale and the lithium abundances, it is also of interest to note 
the effect on an aspect of abundance analysis, specifically ion- 
ization equilibrium often used in the determination of log g. We 
can also make a preliminary investigation into constraints we can 

place on the value of Sh from our results. 

It has been noted previously dGehren et al.l 12001) that Fen 
lines are relatively unaffected by NLTE. In this work we have 
also found this to be the case with values for Fen abundance 
corrections of order 0.01 dex. Our NLTE calculations induce a 
mean difference between AA(Fei) and AA(Fen) of 0.39 dex for 
Sh = and 0.27 dex for Sh = 1 • Knowing that a 0. 1 dex change 
in log g induces a difference of 0.05 dex between Fe i and Fe n 
abundance, for there to be ionization balance, one would need 
a correction of ~ + 0.8 dex and + 0.5 dex in log g for Sh = 
and 1 respectively. That is, due to overionization, forcing ion- 
ization balance for metal-poor dwarfs under LTE calculations 
would give log g values too low by 0.8 dex (Sh = 0) or 0.5 
dex (Sh = 1)- LTE calculations for HD 140283 have occasion- 
ally y ielded gravities lo wer than the hipparcos gravity by ~ 0.3 
(e.g. lRvan et aljfl996bl) . and for a selection of 13 halo main se- 
quance turnoff stars with hipparcos parallaxes Ryan et al. (2009 - 



in preperation) determine a mean difference of 0.2 in log g com- 
pared to LTE ionization balance. These differences are less than 
what we compute for Sh = 1, and suggest that for the model 
atom we are using, the choice of Sh = 1 may underestimate the 
role of collisions with neutral hydrogen in diminishing the de- 
partures from LTE for Fe, i.e. that Sh > 1. Whilst w e have not 
attemp ted a detailed derivation of Sh by this method, Korn et all 

(2003) has, arriving at a value of Sh = 3 based on the analysis of 
four halo stars and two others. Our results are broadly consistent 
with their conclusion. 

7. Conclusions 

We have discussed the processes of NLTE line formation of Fe 
lines. Here, we have shown the challenges posed by such calcu- 
lations and the uncertainties that still arise, in particular due to 
the unknown magnitude of H collisions. As there is at present 
no better theoretical or experimental description of the role of H 
collisions, one obvious next step would be to tie down the value 
of Sh for metal-poor stars, for example by forcing the equality 
of hipparcos gravities an d those determine d by ionization equi- 
librium by changing Sh dKorn et alJl2003h . For this reason we 
have discussed the effect of NLTE corrections on the ionization 
equilibrium and the magnitude of the effect on log g. 

Six of the original program stars from Paper I have been 
analysed to calculate the effects of NLTE on the T e $ scale de- 
rived from Fe i lines via excitation equilibrium. We have found 
that the effect of the correction is to cause an increase in T e g 
ranging from 2 K to 150 K for S H = and 41 K to 122 K for 
Sh = 1. There is one exception; the star CD-33°1 173 has a neg- 
ative correction (-93 K) for the Sh = case. This may be due 
to the limited number of Fe lines available for this star, but also 
emphasises the intricacies of NLTE work which make it difficult 
to make reliable generalisations. 

Our new temperatures have been compa red to the pho- 
tometric temperatures of [ Rvan et alJ (Il999l) . the IRFM of 
Melendez & Ramirez (2004), and the Balmer line wing method 
oflAsplund et al. ( 2006). We find that the NLTE temperatures are 
hotter than IRvan et al.l (1 19991) by an average of 132 K for Sh = 

and 162 K f or Sh — 1 ■ Similar r esults are found when com- 
paring against Asp lund et alJ d2006l) with average differences of 
76 K and 110 K for Sh = and 1 respectively. The d ifference 
between our temperatures and the lAsplund et alJ (120061) temper- 
atures may be rem oved if the Balm er line wing method suffers 
from NLTE effects (Barklem 120071) . or the effects of granulation 
are properly described. We find that even with N LTE corrections 
we are unable to match the high r e ff's of iMelendez & Ramirez! 

(2004) . However, it has been acknowledged that their temper- 
atures suffer from systematic errors (Melendez 2009 - private 
communication) and a revision of their temperature scale is un- 
der way. 

With our new T e s scale, new Li abundances have been cal- 
culat ed. This has led to an increase of the mean Li abundance 
from lHosford et a il (120091) to values of 2. 19 dex with a scatter of 
0.07 dex and 2.21 dex with a scatter of 0.06 dex for Sh = and 

1 respectively, both of which lie well below the val ue of 2.72 
dex inferred from WMAP+BBN (ICvburt et alJ l2008). This has 
shown that systematic errors in the T e g scale of metal-poor stars 
are not the cause for the discrepency. 
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Fig. 1. Departure coefficients for all lower levels (left) and all upper levels (right) of the lines we have studied in the star HD140283. 
A:S H =0,B: S H =0.001,C: S H = 1 
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Fig. 2. Source function,^;, mean intensity,/,,, and Planck function, B v , for two lines whose lower level is close to the ground state 
- X = 0.96 eV (left panel) and 0.12 eV (right panel), for Sh = 1, labelled with the lines characteristic formation depth log T5000 at 
which r v = 2/3. 



